Observation of Nonspreading Wave Packets in an Imaginary Potential 



R. StiitzleQ M.C. Gobel, Th. Horner, E. Kierig, I. Mourachko, and M.K. Oberthaler 

Kirchhoff-Institut fiir Physik, Universitdt Heidelberg, 
Im Neuenheimer Feld 227, D-69120 Heidelberg, German^ 

M.A. Efremov,! M.V. Fedorov,i V.P. Yakovlev,^ K.A.H. van Leeuwen,^ and W.P. Schleich'' 

^General Physics Institute, Russian Academy of Sciences, 38 Vavilov Street, Moscow, 119991 Russia 
^Moscow Engineering Physics Institute (State University), 31 Kashirskoe shosse, Moscow, 115409 Russia 
^Emdhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands and 
*^Abteilung fiir Quantenphysik, Universitdt Ulm, D-89069 Ulm, Germany 
(Dated: February 1, 2008) 

We propose and experimentally demonstrate a method to prepare a nonspreading atomic wave 
packet. Our technique relies on a spatially modulated absorption constantly chiseling away from an 
initially broad de Broglie wave. The resulting contraction is balanced by dispersion due to Heisen- 
berg's uncertainty principle. This quantum evolution results in the formation of a nonspreading 
wave packet of Gaussian form with a spatially quadratic phase. Experimentally, we confirm these 
predictions by observing the evolution of the momentum distribution. Moreover, by employing 
interferometric techniques, we measure the predicted quadratic phase across the wave packet. Non- 
spreading wave packets of this kind also exist in two space dimensions and we can control their 
amplitude and phase using optical elements. 
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Nonspreading wave packets have attracted interest 
since the early days of quantum mechanics. Already in 
1926 Schrodinger |1| found that the displaced Gaussian 
ground state of a harmonic oscillator experiences con- 
formal evolution because a classical force prevents the 
wave packet from spreading. Even in free space the 
correlations between position and momentum stored in 
an initially Airy-function-shaped wave packet can pre- 
vent spreading 2J. Here we propose and experimentally 
observe the formation and propagation of nondispersive 
atomic wave packets in an imaginary (absorptive) poten- 
tial accessible in atom optics ^, ii4„ . Although there 
is no classical force, there are correlations continuously 
imposed by Heisenberg's uncertainty relation resulting in 
the stabilization of the wave packet. 

Localized wave packets due to stabilization are well 
known in the context of periodically driven quantum sys- 
tems (6j and studied with increasing interest for electronic 
wave packets in Rydberg atoms 0, S S 0|- Our ap- 
proach to create nondispersive atomic wave packets re- 
lies on three ingredients: (i) an absorption process 0| 
cuts away the unwanted parts of a broad wave creat- 
ing a packet that is continuously contracting in position 
space, (ii) this process leads due to Heisenberg's uncer- 
tainty relation to a broadening in momentum space and 
consequently to a faster spreading in real space, and (iii) 
the absorptive narrowing and the quantum spreading are 
balanced, leading to a nonspreading wave packet. In the 
following we will refer to such a wave packet as Michelan- 
gelo packet 0. 

Complex potentials for matter waves [l^ emerge from 
the interaction of near resonant light with an open two- 
level system shown in Fig.Q^a). For a standing light wave 



tuned exactly on resonance an array of purely imaginary 
harmonic potentials arises. When the Rabi frequency 
is of the order of the excited state linewidth F the local 
saturation parameter |17o sin(/ca;)/F|, and thus the up- 
per level population, is of the order of unity except in a 
small vicinity of the field nodes. Consequently, our sys- 
tem decays approximately with the rate F. Therefore, 
in the time domain t ^ \/T the atomic wave function 
vanishes almost everywhere, except in small vicinity 5x 
of the field nodes. Here the saturation is small and our 
open system decays with the rate {Vl^kSxY /T <C F. We 
estimate the time dependent size 5x from the relation 
(nok5xft/T - 1 and find 5x{t) - (T /t)^/^ / {k^o). 

The decrease of 5x{t) is accompanied by an increase of 
the width Spit) ~ l/5x{t) in momentum space |l4j lead- 
ing to spatial spreading of the wave packet. Because of 
competition of the two processes - absorptive contraction 
and quantum spreading - the width 5x{t) reaches its min- 
imal stationary value 5xq. In this asymptotic regime the 
rate dx(t)/t of absorptive contraction is obviously bal- 
anced by the rate 6p(t)/M of quantum spreading which 
yields the characteristic time to = I/clJo ^o^O^/^rY^^ 
and the stationary width Sxq ^ (Mcjo)^^/^ with the re- 
coil frequency Ur = k'^/ (2M) <C F. 

The experiments are performed with a slow atomic 
beam of metastable argon {v = 50 m/s) produced with 
a standard Zeeman slower. The brilliance of the beam is 
significantly enhanced with a 2D-M0T setup J^j. The 
final coUimation necessary for coherent illumination is 
obtained by two slits (25 /im and 10 /im) within a dis- 
tance of 25 cm. Applying a Stern-Gerlach magnetic field 
we select the atoms in the internal state Iss (J = 2, 
nij — 0). The imaginary potential is realized with a 
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FIG. 1: Formation of a nonspreading Michelangelo wave 
packet for the center-of-mass motion of an open two-level 
atom (a). The resonant interaction with a standing light wave 
(b) leads to an array of harmonic imaginary potentials. The 
normalized diffraction efficiencies derived from the momen- 
tum distributions (inset) approach a steady state as a func- 
tion of the interaction length Az demonstrating the successful 
realization of stationary wave packets. The solid curves result 
from a numerical integration of the Schrodinger equation [T3| 
with the Rabi frequency Qq = 0.4r. The dashed lines cor- 
respond to the Raman-Nath approximation, revealing that 
the interplay between absorption and quantum spreading is 
essential for obtaining a steady state. 



circularly polarized standing light wave by retroreflect- 
ing a laser beam resonant with the Is5-2p8 transition 
(801 nm). This setup realizes to a very good approxi- 
mation an open two-level system since only 16% of the 
excited atoms fall back to the initial state (in contrast to 
32% without magnetic state selection). In order to con- 
trol the interaction length Az, the laser beam passes an 
adjustable slit. By imaging the slit onto the retroreflect- 
ing mirror we avoid the spoiling effect of light diffraction. 
The detection of the metastable argon atoms is achieved 
by a microchannel plate detector allowing for spatially 
resolved single atom detection utilizing their internal en- 
ergy (12 eV). Since the transverse coherence length of 
the incoming atomic beam is much larger than the opti- 
cal wave length, the outgoing wave function is a coherent 
array of single Michelangelo wave packets, resulting in 
constructive interference in certain directions. The spa- 
tial resolution ~ 50 fim of our atomic detector and the 
free flight distance ^ 0.5 m guarantee clearly resolving 
the resulting atomic diffraction pattern in the far field. 

The diffraction efficiency is deduced by summing up 
the detected number of atoms in angular windows as 
indicated in the right inset of Fig. ^ After their ini- 
tial dynamics the wave packets, i.e., the diffraction ef- 
ficiencies do not change giving evidence to the forma- 
tion of Michelangelo wave packets Our numerical 
simulations (solid line) of the open two-level Schrodinger 



equation take into account the longitudinal as well as 
the transverse velocity distributions Avi = 10 m/s and 
Avt ~ 7 mm/s of the experiment. For fig = 0.4r we 
have a very good agreement with our experimental find- 
ings. Since this agreement depends critically on the Rabi 
frequency we can determine its absolute value. It is con- 
sistent within a factor of 2 both with a rough estimate, 
using the power measurement of the incoming light beam, 
and with the overall absorption of the atomic beam. 

In order to stress that the interplay between absorptive 
narrowing and the quantum spreading is crucial for the 
formation of the Michelangelo packet, we have included 
the result of the Raman-Nath approximation (dashed 
lines). Since this approach is only valid as long as quan- 
tum spreading is negligible, it fails to predict the resulting 
dynamics after the characteristic time to- 

According to the arguments given above, Michelan- 
gelo wave packets emerge after a characteristic time 
to ^ 1/f^o- Our experimental results shown in Fig. [21 
confirm the expected scaling with fimm = 0.23r. 

We now show that a Michelangelo wave packet is a 
complex Gaussian wave packet with a quadratic phase. 
For this purpose we recall that the solution of the 
Schrodinger equation 
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cp{x,t) (1) 



for the metastable state wave function ip{x,t) in the 
vicinity of x = 0, where U2{x) — with luq = 

no v/2w.r/r reads 



(f{x,t) = 



/ k/n 
cosh pt 



exp 



ax tanh /3i 
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(2) 



with a = Mloq exp (— i7r/4) and (3 = luq exp (i7r/4). 
Hence, the probability density |(^(a;,i)p is a Gaus- 
sian with the time dependent width Sx{t) = 
[Re{a tanh(/3i)}] ' , which for ujot > 1 reaches its min- 
imal stationary value k6xQ = {ujr^ /^Iq)^^^ . 

In this asymptotic regime Eq. ^ factorizes into a 
product of the time dependent function cosh^^ {(3t), 
showing that the Michelangelo probability density de- 
cays exponentially in time with the rate Fq = uJo/y/2 <C 
F, and the position dependent complex Gaussian 
exp(— aa;2/2) which contains the quadratic phase 0(a;) = 
MuJox"^ /\/8. A Fourier transformation of this wave 
packet with the stationary width 6x0 , yields the asymp- 
totic behavior of the diffraction efficiencies shown in 
Fig- n by the dashed-dotted fines and is in perfect agree- 
ment with our experimental findings. 

The predicted phase 4>{x) of the Michelangelo packet 
can be deduced from the phases of the observed diffrac- 
tion orders where the phase of the nth order with respect 
to the zeroth order is (j}{n) = — 2(wrF/fiQ)^/2n2 = 02n^- 
To measure the relative phases we realize a compact 
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FIG. 2: Experimental verification of the scaling law to ~ 
1/Qo connecting the characteristic time to = zo/v when 
Michelangelo wave packets form and the Rabi frequency Qo- 
The line is a guide to the eye. We measure the zeroth order 
diffraction efSciency as a function of Az (inset) for different 
Rabi frequencies. The crossing point between the linear ex- 
trapolation of the short and long-time limits yields zq. 

interferometer setup shown in Fig. A thin near- 

resonant probing standing light wave (waist 30 fim) is 
placed directly behind the array of harmonic imaginary 
potentials. The wave function amplitude in each output 
direction is given as a superposition of different diffrac- 
tion orders of the Michelangelo packet. By changing 
the relative phase between the two standing light waves 
we can measure an interference pattern and thus de- 
duce the phase evolution as a function of the interaction 
length Az. 

The interferometric setup employs a probing standing 
wave at 801 nm realized by beams impinging on the mir- 
ror under an angle of 10°. Thus, moving the mirror allows 
us to scan the relative phase between the probing and 
the absorptive light wave (beating period 25 /im). The 
presence of a magnetic field in the interaction region en- 
ables us to realize a detuned (8 MHz) probing wave using 
the same laser for both standing light waves but differ- 
ent circular polarizations. By detuning the probing light 
wave the total flux through the interferometric setup is 
significantly increased in comparison to an exactly reso- 
nant probing field. 

In order to deduce the absolute value of 02 we eval- 
uate the interferometer output in the direction of the 
third diffraction order. For our experimental parameters 
this beam is always a two-beam interference of the first 
and second diffraction order of the array of Michelangelo 
packets. In contrast, the output in lower diffraction order 
directions is the result of multiple-beam interference and 
does not allow us easily to deduce the involved phases. 

In order to find the phase difference 0(2) — 0(1) we 
have to eliminate the offset phase arising mainly from 
the fact that the probing light field is not infinitely thin. 
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FIG. 3: Measurement of the phase of a Michelangelo wave 
packet using an interferometric setup (a) consisting of the ab- 
sorptive and probing standing waves. The inset in (b) shows 
typical interference patterns for different output directions ob- 
tained by scanning the relative position of the second (thin) 
standing light wave for a given interaction length Az. For 
large values of Az the phase shifts (b) of the different inter- 
ferometer outputs relative to the zeroth order level off, indi- 
cating stationary phases of the wave packet. 



For this purpose we take the difference between the mea- 
sured phase in the long-time limit of the absorptive wave 
(Az > 400 /im) and the phase for the experimentally 
achievable shortest interaction length (50 /im). For the 
Rabi frequency flo = (0.23 ± 0.02)r we find the exper- 
imental value 10(2) — 0(1)| — 1.70 ± 0.17, which is in 
agreement with the prediction of the numerical integra- 
tion 0(2) - 0(1) = 302, that is |02| = 0.57 ± 0.1. More- 
over, the characteristic length zq ^ 400 /im for leveling 
off the phases coincides with the one for leveling off the 
diffraction efficiencies. Furthermore, by increasing the 
Rabi frequency to ^Iq = (0.4 ± 0.05)r we experimentally 
deduce |02| = 0.32 ± 0.08, which is in very good agree- 
ment with the prediction of the numerical integration 
|02| = 0.27 ±0.04. 

So far we have concentrated on wave packets in D = 1 
spatial dimensions. A straightforward generalization to 
D — 2 relies on two orthorgonal linear polarized standing 
waves interacting with the appropriate atomic transitions 
and leads to the potential —iM{uj1x^ -\-L0yy'^)/2 near the 
nodes. The frequencies lOx and Uy depend on the field in- 
tensities. A nonorthorgonal configuration provides even 
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FIG. 4: Probability density |(/j(a;)]'^ (solid line) and abso- 
lute value of phase |<?f>(2;)| (dotted line) of Michelangelo wave 
packets in the Potential U2n- 



additional parameters to control the form of the emerging 
two-dimensional Michelangelo wave packet. 

We emphasize that Michelangelo wave packets are not 
restricted to the Gaussian form, Eq. originating from 
the quadratic potential U2 in Eq. Indeed, with an 
appropriate mask 18] we can create almost any behavior 
of the mode function close to the node, leading, for ex- 
ample, to a power law potential U2n{x) = {^lQ/T){qx)'^''^ . 
Here q <^ k determines the characteristic width of U2n- 

The Michelangelo wave packets shown in Fig. ^ for 
n=l,3, and 5 are the "ground" state eigenfunctions of 
the corresponding stationary non-Hermitian Hamiltoni- 
ans and can be obtained numerically. In the asymptotic 
regime only these functions survive because their com- 
plex energy "eigenvalues" have the smallest imaginary 
parts. Moreover, applying the general arguments above 
to the case of U2n yields the following characteristic time 
and width: 



1 



2/(n+l) 



r 



l/(n+l) 



and qSxQ 



l/2r, 



(3) 

where cJ^ — <Z^/(2Af). These scaling behaviors have been 
confirmed by numerical integration of the corresponding 
Schrodinger equation. We note that for n = 1 these 
expressions reduce to the ones of the previous case. For 
n 00 the potential U2n takes on the shape of a box, to 
is independent of flo, and Sxq is solely given by q. 

In conclusion we present a new class of nonspread- 
ing wave packets resulting from the interplay between 
absorptive narrowing and quantum spreading. The de- 
veloped theoretical description explains the experimental 
observation of both the phase and the amplitude of the 
wave packet quantitatively. The experimental realization 
of imaginary potentials strongly relies on spontaneous 
decay processes. Nevertheless, we show that coherence 
is maintained and can even be employed for deducing 
the phase of the Michelangelo packets. Since the wave 
packet arising in the long-time limit is weakly dependent 
on the initial wave function, this process is a robust tool 
for generating wave packets with well-defined amplitude 
and phase for further experiments. 
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